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Abstract 

In this paper we present calculations on the electronic band struc- 
ture of a two-dimensional lateral superlattice subject to a perpendicular 
magnetic field by employing a projection operator technique based on the 
ray-group of magnetotranslation operators. We construct a new basis 
of appropriately symmetrized Bloch-like wavefunctions as linear combi- 
nation of well-localized magnetic- Wannier functions. The magnetic field 
was consistently included in the Wannier functions defined in terms of 
free-electron eigenfunctions in the presence of external magnetic field in 
the symmetric gauge. Using the above basis, we calculate the magnetic 
energy spectrum of electrons in a lateral superlattice with bi-directional 
weak electrostatic modulation. Both a square lattice and a triangular one 
are considered as special cases. Our approach based on group theory han- 
dles the cases of integer and rational magnetic fluxes in a uniform way 
and the provided basis could be convenient for further both analytic and 
numerical calculations. 
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1 Introduction 



Lateral semiconductor superlattices are realized by modern lithographic tech- 
niques to impose a periodic pattern onto a two-dimensional electron system 
(2DES) at a semiconductor heterostructuro. This pattern can be a uni- or bi- 
directional periodic potential modulation of different strength allowing for a 
variety of artificial periodic structures. [1, 2] Magnetotransport measurements in 
such structures have revealed novel oscillations in the conductivity, the Weiss 
oscillations which are superimposed on top of the well-known Shubnikov-de 
Haas oscillations. [3] Besides Weiss oscillations (a commensurability effect), novel 
magneto-resistance oscillations with 1/_B— periodicity have been observed in 
short-period superlattices due to magnetic miniband energy spectrum. [4] The 
above developments have renewed the interest in the problem of Bloch electrons 
in a magnetic field, which was extensively studied in the past by renowned 
theorists [5] but it is experimentally accessible now with these artificial periodic 
structures. [6] 

The first theoretical approach to the problem is due to Peierls[5] and it is 
based on the formulation of an effective single-band Hamiltonian, arising from a 
tight-binding dispersion relation E{k), through the substitution E{k + eA/h), 
where A the vector potential. When combined to a semiclassical picture [7], 
the effective Hamiltonian theory leads to the idea of quantized magnetic orbits, 
which can be commensurate with the lattice period. Within the framework of 
effective single-band Hamiltonian theory, further investigations lead to butterfly- 
like patterns for the energy spectrum known as Hofstadtcr spectrum. [8, 9] This 
spectrum reflects, as a consequence of Harper's equation[10], the splitting of the 
single Bloch-band into magnetic subbands, according to the number of mag- 
netic flux quanta piercing the unit c;ell of the lattic;c. Some recent experimental 
results[ll] supports the existence of the Hofstadter's spectrum in high magnetic 
field regime. Furthermore, an experimental realization of the Hofstadter's spec- 
trum was achieved considering the transmission of microwaves [12] and acoustic 
waves [13] through an array of scatterers. 

In recent years special attention has been paid to investigate the energy 
miniband structure[14]as well as the transport properties of two-dimensional 
short-period lateral supperlattices in the presence of a perpendicular magnetic 
field B in the case of weak modulation. [15, 16] Some calculations of super- 
lattice electronic structure under the above conditions have used perturbation 
theory[17, 18] or a tight-binding scheme[19] and look for solutions of the one- 
electron Hamiltonian that are linear combination of the unperturbed extended 
Landau wavefunctions. Another approach to the problem is based on Ferrari's 
wavefunctions[20], that is, wavefunctions constructed out of coherent states as- 
sociated to a von Neumann-type lattice. [21] Besides the effect of electric mod- 
ulation, the effect of magnetic modulation on the superlatticc energy spectrum 
has been studied[22, 23] and a formal connection between these two has been 
investigated. [24] 

In this paper we present calculations on the electronic band structure of 
a lateral superlattice subject to a perpendicular magnetic field by employing 
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a projection operator technique in the context of magnetic translation group 
(M.T.G) theory. We use a new set of magnetic-consistent Wannier functions 
given in terms of Schrauben functions introduced by Jannussis[25] in order to 
construct Bloch-hkc cigcnfunctions that reflect the crystal symmetry of the su- 
perlattice. Schrauben functions describe free electrons in a uniform magnetic 
field in the symmetric gauge and have been used in the past for magnetocon- 
ductance calciflations.[26, 27] 

The paper has the following structure: In Section 2 the model is formulated 
and the projection operator of the M.T.G corresponding to the problem is intro- 
duced. In Section 3 we construct a basis of appropriately symmetric Bloch-like 
wavefunctions as linear combination of Schrauben functions. In Section 4, the 
electronic spectrum of the 2DES, under weak periodic modulation, is obtained. 
In Section 5, we present specific results for periodic modulations with rectan- 
gular or triangular symmetry. The density of states is calculated as well. We 
conclude with a summarizing Section 6. 



2 Projection operator of magnetic translation 
group 

To describe the electrons in the conduction band of a lateral superlattice at the 
AlGaAs-GaAs interface in a constant perpendicular magnetic field B = Bz, 
we employ a model of strictly two-dimensional electron gas (2DEG) with a 2D 
periodic potential modulation V{r), r = {x,y). 

Following the common practice, we adopt an effective mass approximation, 
to take into account the effect of the crystalline atomic structure over the charge 
carriers. The system is described by the well-known single-particle Hamiltonian 
of Bloch electrons in a magnetic field: 

W=^(P + '=-4)Vy(r) (1) 

where we neglect the Zeeman splitting and spin-orbit interactions which are 
very small for GaAs systems. 

The Hamiltonian (1) does not commute with the usual translations operators 

f{R)=exp{-^R-p) (2) 

where R is any translation vector, but commutes with unitary operators of the 
form 



fm{R)= exp 



(3) 



introduced by Brown[28] and Zak[29], if the vector potential A fulfills the con- 
dition 

1^ + 1^=0 for j,fc = 1,2,3 (4) 
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This relation holds, for the symmetric gauge A{r) ~ ^{B x r), which was 
used by both authors. It is worth noting that introducing local gauge one can 
consider any vector potential. [30, 31] 

The above operators (3) form the magnetic translation group (M.T.G). This 
is in fact a ray group[32], because the product of two elements of the M.T.G 
yields another element, multiplied by a constant phase 



Tra{Rl) Tm{R2) = Tm{Rl + -R2) exp 



--R,-A{R^) 



(5) 



where i?i and i?2 are translation vectors. As a consequence of (5), we can verify 
the commutation relation 



(i2l), Tm{R2) 



= 2i sin 



--B.{R,xR,) 



Tm{Rl+R2) (6) 



where the symmetric gauge for the vector potential has been used. 

Therefore, the magnetic translation operators commute with the Hamilto- 
nian, but in general not with each other. Only in the case that the phase 
constant in (6) is a multiple of tt. 



2h 



B- {Rix R2) =Trq, g € N 



(7) 



we get fm{Ri), TmiRi) = and the operators form an Abelian group. The 
condition (7) is compatible with the demand that the eigenfunctions of the 
system satisfy the Born-von Karman periodic boundary conditions, e.g. are 
periodic under magnetic translations corresponding to the full lattice size. 

The term B ■ x i?2) in (7) is the magnetic flux through the cell spanned 
by lattice vectors Hi and R2 and can be written as an integer multiple {p) of 
the flux per unit cell $ = i>o-B • (oi x 02) , where $0 = is the magnetic flux 
quantum. 

Consequently, the condition(7) is fulfilled each time the number of magnetic 
flux quanta through the unit cell is a rational number: 



$0 



_9 
P 



with q, p mutually prime integers 



(8) 



thus, imposing a commensurability relation between the magnetic field intensity 
and the lattice size. 

After a brief preamble on the M.T.G, we follow Brown's formulation [29] 
to introduce the projection operators of M.T.G. We consider the irreducible 
representations of the M.T.G labelled by the magnetic crystal momentum vector 
q, restricted to a single magnetic Brillouin Zone (MBZ), or to 1/px 1/p part of 
it if p ^ 1. The q — th irreducible representation is given by a familiar relation 



D'^{R) = exp ( - • R) D°{R) 



(9) 



where R = n\a\ + n-ia-i, with a\, a-2 the primitive lattice vectors and q = 
qibi + 92^2 with bi, 62 are the primitive reciprocal-lattice vectors. The vectors 
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ai,2 and bi,2 are related in the usual fashion Oj • bj = 2'K5ij and ni,2, qi,2 are 
integer indices. 

The p X p matrix D^{R) for any translation R is generated from those 
corresponding to the primitive lattice vectors 



-C'fefe'(ai) = ^fefc' exp 



2ni(k~l)- 
P. 



cmodp 



(10) 



and the ray-group multiplication law (5) yields to 



rrnodp 
k^k' —n-2 



exp 



m—ni {n2 + 2(A; — 1)^ 



with 



0,1,2,. 



■p-l 



(11) 

where the modulo Kronecker delta, 6™°'^^ equals 1 if i and j are equal or differ 

by a multiple of p and vanishes otherwise. 

Using (9) and (11) we obtain the following form for the irreducible represen- 
tations 



cmodp 

^ ^k,k'-n2 



2Tri(niqi + n2q2^ 



exp 



-ni(n2 + 2(fc- 1)) 
(12) 

We can see that increasing any of q\, q2 by 1/p produces an equivalent represen- 
tation. Thus the values of q have to be restricted to a single magnetic Brillouin 
Zone: < qi,q2 < ^/p- The rank of matrices implies the existence of j3-partner 
functions transforming according to different rows of the same representation 

^kk' 

We introduce the projection operator of M.T.G. of Bloch electrons in a 
uniform magnetic field[34, 35] 



P^k'iR)=E[^tk'm*Tm{R) 



(13) 



R 



which projects out the component belonging to fc'-th row of the representation 
q and then generates the partner belonging to fc-th row. Using the expression 
for the magnetic translation operator (3), Eq. (13) leads to 



Pkk' (R) = H exp (iq-R) exp 



R 



R (p- 



2h 



BxR) 



where 



Cr =6 



cmodp 



k,k' —n^ 



exp 



— iir-ni (ji2 + 2{k — 1) j 



(14) 



(15) 



We finally remark that acting the projection operator P'^{R) to an arbitrary 
function we can find a partner function of the q-th representation of M.T.G. In 
the next section, we find a basis of Bloch electrons in a magnetic field by acting 
the found projection operator (14) onto a basis of free electrons in a magnetic 
field. 
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3 Localized wave- functions of Bloch electrons in 
a magnetic field 

Landau functions, eigenfunctions of angular momentum and localized functions 
based on coherent states [21. 36] arc usually used to describe free electrons mov- 
ing in a plane under the influence of a magnetic field. The coherent states[37] on 
a von Neumann lattice are the most localized eigenfunctions for an electron in 
magnetic field. For such a lattice, with a single fiux quantum per unit cell, the 
area of the unit cell equals 2nP, where I = {h/eBy^^ is the magnetic length. 

Jannussis[25] working in the symmetric gauge, introduced a different basis 
of wave functions for free electrons in uniform magnetic field, who called them 
Schrauben functions: 



(4 + + ik-r 



(16) 



where n e N is the Landau index, k the electron wavevector and Kx = kx+y/2P, 
Ky = ky — x/2P, The normalization coefficient is given by 



j_{2py>_ 

2wP n! 



(17) 



The above functions have several important properties, some of which are: 

i)They form a basis in the Hilbert space of free electrons in a uniform mag- 
netic field. 

T.rnA^)Mr') = S{r-r') (18) 

n,k 

(19) 



(20) 



ii)They are fully localized around the centers 
and in the two-dimensional plane take the symmetric form : 



i^n,k{'^) =C„ 



{x - Xmf - i{y - y-, 
1 



• exp 



4/2 



(X-Xm) +{y-ym) +^iXmy-Xym) 



(21) 



Also, the probability density p = ip*ip has a rotational symmetry about the 
centers {Xm,ym) and takes the form of the Poisson distribution 



1 jrll^l' 
2ttP n\ 



■exp - 



2f 



(22) 



where = (a; - Xmf + {y - Vmf 
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iii) They obey Harper's symmetric condition[10], i.e. for a lattice translation 

R 



i^n,k{r + R)= exp 



A{R) r + ik r 



K = k--A(R) (23) 

h 



It has been shown[38] that a Wannier function ui„(r — r^) for electrons in 
a magnetic field can be obtained directly form the Schrauben function if the 
wave- vector k is substituted by the vector potential —eA{rm)/fi 

In what follows, we construct a basis for two-dimensional Bloch electrons 
in a uniform magnetic field applying the projection operator method described 
in Section 2. In order to find the principal partner function we choose the 
projection operator P'}^{R) (see Eq. (14)), projecting onto the first row of q-ih 
representation. 

Acting the projection operator Eq. (14) onto the functions V'ti,o('") of the 
basis (16) we find the following partner functions of the q-th representation of 
the M.T.G. 



(r) 



R 



- ^ Cr exp [iq ■ i?)w„(r - R), ko = -^A{R) 



R 



(24) 



where Wn{r — R) are the magnetic Wannier functions of free electrons in a 
uniform magnetic field [38], 



Wn{r - R) 



:^[(a;-niai)+i(y- 72202)] 
1 



• exp 



4/2 



(x-niai) +{y-n2a2y 



+ 



:{yniai - 0:71202) 
(25) 



where ai, 02 are the 2D-lattice periods along x and y directions respectively 
and the normalization constant C„ is given by Eq.(17). 

Therefore, Eq. (24) represents modified Bloch functions which are formally 
expanded in the corresponding magnetic Wannier functions, that is, functions 
strongly localized on the lattice sites. Note that, for commensurate magnetic 
fields (7) one can recover Bloch's theorem 



f{a)^n,qif) = exp ( - • a)*„,g(r) 



(26) 



with the values of the magnetic crystal momentum q = {qx,<ly) restricted to 
the first magnetic Brillouin zone. Also, ^'„,q with diflferent values of q are 
orthogonal since they belong to different irreducible representations of M.T.G. 
Finally, after some straightforward steps, we arrive at the following expression 
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for the eigenfunctions Eq. (24) 



^n,q{r) = ^ Cnin2Cn exp 



i{qxniai + qyn2a2) 



2p[{^ - niai) + i{y - 71202) 



• exp 



4/2 



{x - niaif + {y- 71202)' 



2/2 



(ynifli - xn2a2) 



(27) 



where C„i„2 = exp {—mnin2q/p) with n2 integer multiple of p. The values 
of crystal momentum q = {qx,Qy) are determined from the periodic boundary 
conditions for ^'„^q(r) with respect to magnetic translations with periods ai 
and a2. 

Figure 1 shows the contour plots of the modulus |'J'o,q(T')| for the case of a 

square periodic lattice (ai = 02 = lOOnm) in the presence of a perpendicular 
magnetic field that corresponds to two flux quanta per unit cell, for different 
values of the crystal momentum q in the first Brillouin zone. It is seen that 
the wavcfunctions possess the full symmetry of the lattice and have regularly 
distributed zeros as it is expected from the translational invariance of the system. 



4 Electron energy spectrum under weak peri- 
odic modulation potential 

Now we can find the energy spectrum of Bloch electrons in a constant magnetic 
field and weak periodic modulation potential V{r) = V{r + R) using the basis 
||n, <7)} introduced in the preceding Section 3. 

It is well-known [25] that the matrix elements of the first part of the Hamil- 
tonian (1) which describes free electrons in a uniform magnetic field are given 

by 



n,q 



1 / 



2m* 



1 



- (p + eA) n', qj = {n+- )hujc 5nn' (28) 

where lo^ — eB/m* is the cyclotron frequency and n labels the Landau levels. 

In order to find the corresponding energy spectrum of Bloch electrons we first 
evaluate the matrix elements of the periodic potential. Since V{r) is periodic 
with respect to R it can be expressed in terms of a Fourier series 

V{r) = ^Vg exp [iG ■ r) (29) 

G 

with the Fourier coeSicients given by 



Vg = \ I d'^rV{r)exp{-iG-r) 

J A 



(30) 



where A stands for the area of the unit cell in the two-dimensional lattice and 
G is a reciprocal lattice vector. 
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In the two dimensional lattice {x, y) the potential of Eq. (29) is written as 

'27r 2tt 



V{x,y) ^ yVg^^g^exp i[ — gix-\ g2y 



31 S2 



(31) 



with 5i, 32 integers and Oi, a-i the modulation periods in x and y direction, 
respectively. 

Using the basis (27), the matrix elements of the potential F(a;,y) can be 
derived for each crystal momentum q= {qx,(ly)'- 



Vnn' ={n,q 

■ exp 

• exp 

• exp 

where 



V{x,y) n',qj = J2 Qi«2^"i"i^9i92 

9192 nin2 n'^n'^ 

iqxai{ni - n[) + iqya2{n2 - 



("2 - n'2){n2 -n2+ 251^) - ^(ni - n[){ni ~ n[ - 2^2^) 



ITT 



{n[n2 - n'2ni)— + (m + n[)gi + (n2 + 712)92 
3>o 



Inn' {91,92) 



(32) 



ai 



a2 



Inn' {91, 92) ={-^) \ -r \ —92 + ix —91 



a2 



ai 



$0 



(n-n')/2 



■exp{-u/2)L'^7"-'{u), for n>n' 



(33) 



and 



/„n'(5l, 52) =(-!)"+" 



-92 + i\ 



-91 



£0 



(n'-Tl)/2 



■exp(-u/2)L;f-"(u), n'>n 



(34) 



with u = TT{—gi H — -^2)"^ ^^'^ I'mi^) the associate Laguerre polynomials, 
oi 02 ^ 

Furthermore, from the expression (32) one sees that for high magnetic fields 
or if the lattice points are far apart each other, the matrix elements Vnn' ap- 
proach zero unless the following conditions are valid: 



m — Til = 2g2— and n2 — "•2 = "^gi — 



(35) 



Under the above conditions, which mean that the magnetic flux through a unit 
cell of the periodic lattice should be a rational number, the matrix elements 
(32) take the form 



Vnn' = ^ Vg^g^ exp 2i^{92aiqx - 5102%) Inn' {91, 92) 



91 92 



(36) 
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Finally, the first order energy perturbation can be found by calculating the 
diagonal matrix elements from Eqs. (36) and (33); 

Vnn{qx, Qy) = Vg^g^ cxp [2^— (gsaife - 51^2%)] cxp (-u/2)L„(u) (37) 

9192 

Thus, the complete energy spectrum to first order approximation can be written 
as ^ 

En{qx, qy) = {n+- )hu)c + Vnn{qx, qy) (38) 

5 Energy band structure for lateral superlat- 
tices with weak modulation 

In this section, we present analytical results concerning the energy spectrum 
of superlattices with perfect rectangular or triangular symmetry applying the 
formalism of the preceding Section 4. We consider smooth periodic potentials 
described by a few lowest Fourier components with weak modulation strengths. 

It is well-known that the numerical solution of the corresponding tight- 
binding equation[39] for 2D-lattices in perpendicular magnetic field, shows that 
each Landau level splits into several sub-bands when a rational number of flux 
quanta $o pierces the unit cell and that the corresponding gaps are exponentially 
small. In what follows, we assume that the gaps are closed due to disorder in 
samples of not exceptionally high mobilities, that is, we assume that the Landau 
levels are bands neglecting the above fine structure. This assumption is justified 
since the reported experiments [40] did not indicate that this fine structure of the 
energy spectrum was resolved. 

5.1 Lateral superlattices with rectangular symmetry 

We consider a two-dimensional periodic modulation with rectangular symmetry, 
i.e: 

27r 2-77 
V{x, y) = cos {—x) + Vy cos (— y) (39) 

where 14, Vy are the modulation strengths and ai 02 the periods in the corre- 
sponding directions. We assume that the modulation is weak compared to the 
cyclotron energy fkOc so that Eq. (38) is valid and mixing of different Landau 
levels can be neglected. 

In terms of the plane wave expansion Eq. (31), the periodic potential (39) 
can be described by four Fourier coefficients 

V(±i,o) = T4/2 and l/(o,±i) = Vy/2 
Applying Eq. (38) we find the energy spectrum 



E{qx, qy) = {n+l; )hujc + VxF{ux) cos (^^Uxqy)+VyF{uy) cos i^^Uyq^) 

I 'K IT 



(40) 
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where 



F{u) = exp{-u/2)Ln{u), = —{n^o/^) and = — (tt^o/^) 

CLi a2 

Thus, the unperturbed Landau levels broaden into bands, with width equal 
to 2(Vx\Fn{ux)\ + Vy\Fn{uy)\), that oscillates with inverse magnetic flux $o/^ 
. Using the asymptotic form of Lagucrre polynomials in the large n- limit, it 
can be shown that the bandwidth tends to zero when = (A + l/4)7r and 

reaches maximum when 2\/mr = (A — 1/A)w with A an integer. Therefore, the 
energy bandwidths in Eq. (40) oscillate with the Landau index. One important 
consequence of the non-zero bandwidth is finite electron-velocities leading to a 
finite band conductivity[41] which is absent when the modulation is not present. 

In Fig. 2 wc plot the energy spectrum for the first two Landau levels, n = 0, 1, 
modulation strengths = Vy = 0.5 meV and for (a) an integer ($/$o = 2) and 
(b) a rational ($/$o = 5/4) number of flux quanta. The square superlattice 
period has been taken a = 100 nm. It is seen that the energy spectrum reflects 
the symmetry of modulated square lattice. 

The interesting characteristics in the energy spectrum of the weak modulated 
superlattice are also reflected in the density of states (DOS) per unit area deflned 
by 

= ^ E E -^[^ - Eniq.,qy)] (41) 

where = 2 is the electro-spin degeneracy. For the spectrum given by Eq. (40) 
the DOS per unit area takes the following form: 

g{E) =go V / d0{[V,Fn{u,)]^ -[E-{n + 1/2 )huj, - VyF^iuy) cos 0]^}-^/^ 

(42) 

where go = l/(7r^Z^). Figure 3 plots the DOS per unit area for the same values 
of parameters as in Fig. 2. As can be seen, the DOS does not exhibit Van 
Hove singularities at the edges of each Landau band reflecting the 2D nature of 
the electron motion in this band with finite mean velocities. Finally, wc have 
to note that similar results for the band structure of a two dimensional surface 
superlattice have been obtained previously by direct application of perturbation 
theory to a different unperturbed basis of extended Landau wavefunctions. [18, 
39] 



5.2 Lateral superlattices with triangular symmetry 

We now consider a laterally modulated 2DEG with triangular symmetry, i.e: 

27r 27r 47r 

V{x, y) = Vx cos ( — x) cos ( — y) + Vy cos ( — y) (43) 
ai a2 (12 
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where ai = a and a2 = VSa. This periodic modulation can be reaUzed in a 
GaAs/AlGaAs heterostructure by means of a diblock copolymer nanolithogra- 
phy technique. [42] 

The above periodic potential can be described by six Fourier coefficients 

V(i,±i) = V(_i,±i) = and y(o,±2) = Vy/2 

where we assume again the condition of weak modulation Vx,Vy « Huc- Ap- 
plying Eq. (38) we find the energy spectrum 

E{qx, %) ={n + 1/2 )hujc + VxF{ux + Uy) cos {—Uxqy) cos (^^^Wj/^x) 

TT TT 

+VyF{4Uy) COS i^^^Uyq^) (44) 

TT 

where the parameters Ux, Uy and the function Fn{u) are defined in the preceding 
subsection 5.1. 

As in the rectangular case, the Landau levels have broadened into bands with 
a bandwidth equal to 2{Vx + Vy)\Fn{Auy) , that oscillates with inverse magnetic 
flux $o/^ and Landau band index n. In Fig. 4 the energy bands for n = and 
n = 1, are plotted for the same values of parameters as in the rectangular case. 
The energy spectrum reflects the symmetry of the modulated triangular lattice. 

Furthermore, the DOS per unit area of the superlattice with triangular sym- 
metry can be calculated according to Eq. (41) and is given by 

<.27r 



g{E) =go V / d0{[VxF^{ux + %) cos^]^ 



-[E-{n + 1/2 )hujc - VyFn{4uy) cos2(9]^}-^/^ (45) 

where go = l/{Tr^P). The DOS for this case has similar characteristics to the 
one shown in Fig. 3 and is not plotted. 

5.3 Exact energy spectrum 

The approximate analytic energy spectrum of the preceding subsections 5.1 and 
5.2, provides a good insight into the problem and can facilitate the study of 
other properties, e.g., the thermodynamic and the transport properties of the 
superlattice. On the other hand, the exact energy spectrum can be obtained 
by diagonalizing numerically the Hamiltonian matrix produced by means of 
Eq.(36). The diagonalization is carried out by truncating the matrix dimensions 
while including a sufficient number of Landau levels. We were careful to choose 
the rank of the matrix being diagonalized large enough to ensure the convergence 
of our result. In Fig. 5 the exact energy spectrum of the superlattice with square 
symmetry, given by Eq.(40) and the approximate one are plotted as a function 
of the crystal momentum (qx, Qy = 0) for the first ten Landau levels, where the 
modulation strengths are Vx = Vy = 0.5 meV and $/$o = 2. 
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Wc sec that the approximate energy spectrum (soHd curve) agrees well in 
the whole range of crystal momentum with the exact one (dotted line). Similar 
agreement is found for other values of magnetic flux $/$o provided that the 
condition of weak modulation is satisfied. 

6 Concluding Remarks 

In summary, we have studied the energy spectrum of electrons in lateral semi- 
conductor superlattices subject to a perpendicular magnetic field. Such systems, 
in which the simultaneous influence of a bi-directional periodic potential and a 
magnetic field is studied, are realized nowadays in layered semiconductor struc- 
tures with modern technologies of growth and lithography. 

We introduced a projection operator method applied to irreducible repre- 
sentations of the magnetotranslation group in order to construct an appropri- 
ately symmetrized basis of Bloch-type wavefunctions in terms of a new set of 
magnetic-consistent Wannier functions. The well-localized Wannier functions 
in this linear combination, were produced by eigenfunctions of free electrons in 
a magnetic field obtained in the symmetric gauge. 

Using the above Bloch-type wavefunctions we have calculated the magnetic 
energy spcctriim of Bloch-clcctrons in a two dimensional periodic potential ob- 
taining similar results to previous pertubative or numerical approaches. The 
spectrum is expressed as a function of crystal-momenta, defined in a magnetic 
Brillouin zone, which are good quantum numbers of the system. Especially, we 
have presented explicit results concerning the band structure of short-period 
lateral superlattices with both square and triangular symmetry. The obtained 
oscillatory behavior of energy bandwidths, in the weak modulation regime, man- 
ifest themselves in the magnetoresistance measurements where they give rise to 
the so called Weiss oscillations in the dependance of the conductance on the 
magnetic field strength. 

Finally, we note that our approach provides a symmetry-adapted basis- 
functions based on group theory and handles the cases of integer and rational 
magnetic flux in a uniform way. For rational magnetic fluxes it is assiimcd that 
the fine structure of Landau bands can be neglected in samples of not excep- 
tionally high mobilities, since the mini-gaps in a Landau band are closed due to 
disorder. From the presented results, it seems that our basis is appropriate in 
the intermediate regime of magnetic fields where both the lattice period and the 
magnetic length are comparable in magnitude. We also expect that this basis 
is convenient for further both analytic and numerical calculations. 
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Figure 1: Contour plots of the modulus |\E'o,q(r)| for a square lattice with a = 
100 nm and $/$o = 2, for different values of the crystal momentum q in 

the first Brillouin zone: (a) [q^^qy) = (0,0), (b) {qx,Qy) = (""/ajO), (c) 
(&,53/) = (0,7r/a), (d) {qx,qy) = {-k / a,iT / a) 
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Figure 2: Magnetic energy bands n = 0, 1, of the square lattice as a function of 
the crystal momentum (g^, Qy) for (A) $/$o = 2. and (B) for $/>I>o = 5/4.The 
modulation strengths are Vx = Vy = 0.5 meV and the energy eigenvalues are 
plotted in units oieQ = hw. 
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Figure 3: Density of states per unit area of the square lattice for the same values 
of parameters as in Fig. 2. 
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Figure 4: Magnetic energy bands n = 0, l,of the triangular lattice as a function 
of the crystal momentum {qx, Qy) for or (A) $/$o = 2. and (B) for <I>/$o = 5/4. 
The modulation strengths are 14 = ^ = 0.5 meV and the energy eigenvalues 
are plotted in units oi eQ = Tvu). 
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Figure 5: Energy spectrum of the square lattice as a function of the crystal 
momentum {qxT<ly = 0) for the same values of parameters as in Fig. 2. Solid 
lines are for the first order energy perturbation and dotted lines are for the 
exact energy spectrum. Only the lowest ten levels are shown. The energy 
eigenvalues are plotted in units of eo = ^ 
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